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The interplay of dilatonic effects in dilaton cosmology and stochastic quantum space-time defects 
within the framework of string/brane cosmologies is examined. The Boltzmann equation describes 
the physics of thermal dark-matter-relic abundances in the presence of rolling dilatons. These 
dilatons affect the coupling of stringy matter to D-particle defects, which are generic in string 
£SJ \ theory. This coupling leads to an additional source term in the Boltzmann equation. The techniques 

of asymptotic matching and boundary-layer theory, which were recently applied by two of the 
authors (CMB and SS) to a Boltzmann equation, are used here to find the detailed asymptotic 
CN) ' relic abundances for all ranges of the expectation value of the dilaton field. The phenomenological 

q , implications for the search of supersymmetric dark matter in current colliders, such as the LHC, are 

• discussed. 
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I. INTRODUCTION 

To evaluate candidates for cold dark matter (DM) it is necessary to compute relic abundances including physics 
I . beyond the standard model. String-theory considerations provide a natural source of such physics. 

Recently, a global asymptotic analysis was performed on the (Riccati-type) Boltzmann differential equation that 
describes the evolution of the thermal DM relic abundances in an expanding universe [l| . It was shown that boundary- 
layer theory, which makes use of asymptotic matching Q, can give a consistent approximate solution to this Riccati 
equation in two physically interesting cases: (i) standard Friedman- Robertson- Walker (FRW) cosmology Q, and (ii) 
dilatonic string cosmology 0j @. In case (i) the freeze-out and post-freeze-out regions (we emphasize that these 
are regions and not isolated points) for the DM abundances were denned using this novel approach. In case (ii) 
the Boltzman equation of case (i) is modified by the addition of a rolling dilaton source term derivable from string 
theory and proportional to the dilaton cosmic rate The effects of the rolling dilaton on cold DM abundances 
were calculated and it was shown that there is a large-time power-law decay of the DM abundance (with calculable 
■ corrections) . The latter results explain the findings of @ on the dilution of DM relic abundances in the current epoch 
in supersymmetric theories with rolling dilatons. This dilution may significantly affect the available parameter space 
£SJ ■ (after the appropriate cosmological constraints from WMAP 0] are taken into account) and, in turn, may affect the 
t— I | searches for supersymmetry at colliders such as the LHC [1]. 

The analysis cited above does not include the effect of a cosmological background due to effectively point-like defects 
(quantum space-time foam), which are generically found in models based on string theory Dilatons are coupled 
to the foam through the string coupling constant. This foam modifies the effect of the dilaton in the evolution of 
DM and can even dominate asymptotically in the absence of dilaton effects. In our model of space-time foam the 
universe is represented as a brane, with three large spatial longitudinal dimensions, embedded in a higher-dimensional 
bulk space. The "foamy" structures are provided by stringy membrane (D-brane) defects, compactified appropriately 
along extra-dimensional manifolds. From the point of view of a four-dimensional observer the defects appear to be 
point-like (D particles). 1 As the D-brane world moves in the bulk, the D particles cross it and thus appear to the 
four-dimensional observer as stochastic space-time structures, flashing on and off. The stochasticity in target-space 
is attributed to quantum fluctuations of the D particles, viewed as stringy dynamical entities embedded in the bulk 
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space. The dilaton $ directly affects this process because its vacuum expectation value determines the string coupling 
g s . This paper investigates the interplay between the dilaton and a background of space-time defects and their effects 
on the asymptotic behavior of relic abundances. 

In Sec.|TT]we review the main results of [l[ concerning the DM relic density for asymptotically longtimes in standard 
and dilatonic cosmologies. This serves to introduce the powerful technique of asymptotic matching [2j used in Ref. 
In Sec. IIIII we apply the analytical methods of [l| to the case of D-particle stochastic foam (D foam) in the presence 
of relaxing dilatons. The D foam is characterized by a source that differs from the source in dilatonic cosmology. 
Various classes of asymptotic behaviors are determined by the expectation value of <fr. In Sec. IIVI we discuss the 
phenomenology of these models. We address the combined effects of the running-dilaton and D-foam sources on 
the thermal DM relic abundances today and the associated constraints implied by the current LHC phenomenology. 
Finally, a technical discussion of the thermodynamic properties of the various types of universes in the presence of 
sources, which are examined in this article, is given in the Appendix. There, we explain in detail how it is possible 
to define an entropy function that is conserved in the presence of nontrivial source terms in the Boltzmann equation; 
this function allows a thermodynamic interpretation of the respective cosmological equations. 



II. REVIEW OF ASYMPTOTIC SOLUTIONS TO THE BOLTZMANN EQUATIONS FOR RELIC 
ABUNDANCES IN STANDARD AND DILATONIC COSMOLOGIES 

For a DM species X of mass mx the evolution of Y(x) = Af/s, the number density W per entropy density s, is 
governed by the Boltzmann equation [3] 

y'(x) = -Ax— 2 [r 2 (a;)-V c 2 q (x)], (1) 

where x = mx /T is the dimensionless independent variable and T is the temperature. This Riccati equation does 
not include any dilatonic effects of string theory. We are primarily interested in epochs of the universe for which 
mx > T > To, where To is the current temperature of the universe. The integer n — 0, 1, 2, . . . comes from a 
partial-wave analysis of the scattering of DM particles: n = refers to s-wave scattering, n = 1 characterizes p-wave 
scattering, and so on. The parameter A is a dimensionless measure of the scattering of DM particles and is regarded 
as a large number (A > 1). If we parametrize the thermally averaged annihilation cross-section (av) = ao x~ n with 
n = 0, 1, . . . for (s, p, . . . )-wave DM annihilation, and the Hubble parameter as H — H m x~ 2 , then A = <7o "n^c/ H m 
3]. For bosonic remnants the function Y cq (x) is the distribution [ll| 

r°° s 2 

Y cq (x)=A^ ds eVsr ^_ i , (2) 

where A = 0.145g/<7*, g is the degeneracy factor for the DM species, and g* counts the total number of massless 
degrees of freedom Q. 

A closed-form analytical solution to the Riccati equation (fT]) is unavailable, so an approximate heuristic approach 
is customarily used to treat this equation: As the universe cools and x increases, the nature of the solution Y{x) to 
(fT]) changes rapidly in the vicinity of a value x = Xf, the so-called freeze-out point, and as x — > oo the solution Y{x) 
approaches the constant Y^ , called the relic abundance. One approximation is made for x < Xf and another is made 
for x > Xf. The solutions in the two regions are then patched at x = Xt. The value Xt is determined from equating 
the interaction rate of the DM particle and the expansion rate of the universe, a sensible physical criterion. 

This approach gives a reasonably accurate determination of Y^ and, prior to the work of [l|, it has been widely 
adopted [3j. However, this splitting into two regions is only a pragmatic convenience and there is really no precise 
value Xf. Rather, there may be (in a sense to be specified) a freeze-out region. Because the differential equation ([T]) is 
first order, its solution is completely determined by one initial condition, namely Y(0). The usual method of splitting 
|T]) into two approximate first-order equations, which are valid in each of two regions, requires two conditions, an 
initial condition and a patching condition. The value of Xf becomes explicitly involved in the determination of Y^ 
even though the mathematical theory of differential equations does not require this. To avoid this unsatisfactory 
mathematical treatment (which is common in the literature), two of the current authors (CMB and SS) presented in 
[l| a detailed analysis of the associated Riccati equations using applied mathematical methods commonly used in fluid 
mechanics. A key concept is that the freeze-out region can, at least in physically relevant cases, be considered as a 
boundary layer. The solutions in the two regions can then be matched asymptotically. Before reviewing the solution 
of (TT|) for large x we introduce the Boltzmann equation in the presence of a dilaton background. 

In the case of rolling dilaton cosmologies [4( the thermal DM relic abundance is characterized by the presence of 
a linear sink term, which is proportional to the rate of the rolling dilaton field ^ Q. In theories with scale- factor 
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duality [lj . we have 

$(t) = $0 log a(t), (3) 

where a(t) is the scale factor of the expanding universe. In eras where the temperature T satisfies mx > T > To, 

Y'{x) = -Ax-"" 2 [Y 2 (x) - Y c \(x)} + * Y(x)/x. (4) 

Here, $o is a negative dimensionless constant of order 1 that appears in the general expression for the dilaton field 
as a function of cosmic time t. For $o = — <P < 0; the string coupling g s = e* becomes perturbatively small for large 
times and vanishes asymptotically as t — > oo. Thus, the er-model perturbative picture suffices to describe the features 
of cosmology at large-times. As shown in [lj, the presence of the particular dilaton source in (j4]) gives a solution for 
Y(x) whose behavior is qualitatively different from the solution for Y{x) in ([!]). 



A. Boundary-Layer Theory 

Since A is large, the highest derivative in both equations (JTJ and (0]) is multiplied by a small parameter, which 
implies that these equations may be treated by using boundary-layer techniques Q and leads to the concept of a 
freeze-out region as opposed to a freeze-out point [lj]. We rewrite Q as 

jZ'(x) = -x— 2 [x^Z\x) - x+Y&x)] , (5) 

where 

Z(x) = Y(x)x' p . (6) 

The coefficient 1/A of the highest-derivative term is very small. The number of terms on the right side has been 
reduced from three to two; this facilitates asymptotic matching. Outside a boundary layer (the outer region), Z(x) 
varies slowly. Inside a boundary layer, Z{x) varies rapidly. 

We have two outer regions where Z{x) = Z^ x \x) and Z(x) — Z^ 2 \x), respectively. In the left outer region 
Z^(x) w Z cq (x) = x^Y cq (x). To be precise, we write 

oo 

zV{x)~52\-*z£-\z). (7) 

fe=0 

On substituting Z^(x) into (O, we find that 

Z { 1] (x) = Ae~ x x^\ 

Z^ix) = x v+n+2 /2, (8) 
and so on. The entity Xf is defined to be the value of x for which 

zW(x) = z[ 1 \x) (9) 

and is a measure of where the equilibrium region ends. Equation (j9|) implies that 

Xf ~ log(2AA) - (n + 1/2) log (x/). (10) 

This analysis is somewhat simplified (see [l[). The higher order terms in ([7]) are not negligible, but they lead to a series 
with alternating signs that is Borel summable. The Borel sum of the series leads to a multiplicative renormalization 
of A by a factor close to 1. In order to keep the notation simple we have not distinguished A from the renomalized 
A. Solving the equation obtained by replacing in (|10|) the symbol ~ by the equality sign gives Xf. 



Xf = (n + l/2)W 



(2XAy 



-1/2 



1/2 



(11) 



where W(z) is a Lambert function Hence the asymptotic behavior is fully determined in terms of constants 

occurring in the Boltzmann equation. 
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The value Xf lies in the transition region from equilibrium to freeze-out which is interpreted as a boundary layer. 
This interpretation can be justified by the method of asymptotic matching. We define an inner variable X as follows: 



x = Xf + kX. 

Then \X\ can be large compared to 1 but small compared to A. Now, for Z(X) we have 



-Z'(X) = -Xx 



-n-2-<j> 



f 



Z 2 (X)-A 2 x 3 f +2 *e- 2x f 



-Xa 



-n—2 — tf> v 2 



f 



Z\X). 



The exponential term is negligible because Xf ss 25 for typical values A « 10 14 and A w 0.00145. 
From the principle of dominant balance 0] we have 



k = x 



n+2+4) 
f 



/A. 



The solution to ([13)) is 



Z(X) = 1/(X+D), 



where D is a constant of integration. This is the solution in the boundary-layer (or freeze-out) region. 
To the right of this boundary layer there is a second outer region. For large x in this region 



Z'(x) ~ -Xx- ,l - 2 -*Z 2 (x) (x » 1), 



whose solution is 



^post- 



- freeze— out 



(X) 



l/C - Xx-'^-'t'/in + 1 + <t>) ' 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 



where C is an integration constant. 

The behaviors in the equilibrium outer region, the boundary-layer region and the post-freeze-out outer region must 
be asymptotically matched. This matching determines the constants of integration C and D. We first match the 
solution in the equilibrium region to the boundary-layer solution: 



Z 



thermal— c> 



quiUbrium^) ~ 2Ax 3/2+4> e- x - 2^4(2;/ + nXf' 2 ^ e~ x f e~ K 



x 



The factor of 2 is included because two lowest-order terms of the expansion in (|7|) are considered. Noting that k and 
X/xf are small, we get 



^thermal— equilibrium^ 



n+2+4, 
C f 



X(1 + kX) X + Xxj n - 2 -^ 



on using HH). Hence, comparing with (fT5|) . we deduce that 



D = XxZ 



(18) 



(19) 



Similarly, (IT7|) leads to 



from which we deduce that 



ost — freeze— out 



C n+1+0 



Comparing with ([15]) . we get 



^post- 



- freeze— out 



(x)~- 



C {n+l+4)x n f 



X 



(20) 



D C !», i I " 1 " " 



(n + 1 + <j>) x n p 



(21) 
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Finally, from (jT9|) we deduce that 



(n + 1 + 0)x" +2+ ^ 

_ yy / (22) 

A(n + l + <?!> + a;/)' 

The leading behavior for large x in the post-freeze-out region is 

(n + 1 + 0K +2+0 , 

Y ( x ) ~ T7 T-^T • (23) 

w A(n + 1 + </> + £/) v ; 

We denote the solution to ([1} as 1^,5(2), where ns stands for no source. Its asymptotic value for large x is obtained 
from (|2"3"|) by setting = 0. The specific solution for xt in (fTTj) is denoted by Xf ins . 

The above calculation forms the basis of the following analysis that will be given for various parameter ranges and 
sources in the Boltzmann equation. 



III. DM RELIC ABUNDANCES: THE CASE OF A STOCHASTIC STRINGY SPACE-TIME FOAM 

The background of stochastic D-particle foam leads [9j to the inclusion of a positive source T (as opposed to the 
sink in dilaton cosmology) in the standard Boltzmann equation for the thermal relic abundance of the DM species X 
of mass rax ■ In terms of the number density TV it was shown in Ref. |9] that the Boltzmann equation reads 

^+3HAf=T(t)Af+C[f], (24) 

where C[f] denotes the Boltzmann interaction terms and 

a 2 

r(t) = 2Hm x a 4 (t)-^T(9 + 2m x /T) < A 2 >, (25) 
M£ 

where M s is the string mass scale. (The mass of a D-particle defect in the foam is M s /g s Q.) The quantity <C A 2 ^> 
is a dimensionless variable, which expresses the variance in the recoil velocities of the D-particle defects in the foam, 
during their collisions with the DM particles [9j . 

The symbol <C • • • ^> denotes the average over the population of D particles on the three-dimensional-space brane 
world in a given epoch of the universe. The no-force (dust-like) behavior of the D particles, implies the following 
scaling of <C A 2 ^> with the scale factor a(t) of the four-dimensional (brane) universe: 

« A 2 »= <A 2 > a~ 3 (t) = (A 2 ) C - 3 T 3 = (A 2 ) m\ C^ 3 x' 3 . (26) 

Here Co = a(to)7o is a dimensionful constant that appears in the cooling law of the universe; that is, 

a(t)=C Q /T = a(t Q )/(l + z), (27) 

where z is the redshift parameter. The values z = 0, t = t$, and T = To are correspond to the current era. This 
source is positive (in contrast to the sink of dilaton cosmology) and is discussed in a more general framework in the 
Appendix. 

We now discuss the collision term (av) 

number density. Eq. (|24[) now becomes 



(A/" (0) ) -A/" 2 in flU}, where A/" (0) is the equilibrium value of the DM 



Y'(x) = -Ax-™" 2 [Y 2 (x) -Y c 2 (x)] + « A 2 » x 2 (9 + 2x) Y (x) . (28) 

1 m x Mi 



Hence, the Boltzmann equation (|28p becomes 

Y'{x) = -Ax— 2 [Y 2 (x) - Y 2 q (x)] + g 2 (A 2 } (9 + 2x) Y(x)/x. (29) 

s 

There is an implicit dilaton dependence in (|29[) that needs to be made explicit. The string coupling g s is the 
exponential of the dilaton, g s = go exp ((<£)), and so 

gs^goa-^goiCo/rnxy^x-t (30) 
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Hence, for consistency we must incorporate both the dilaton sink and the source induced by D-particle foam in the 
Boltzmann equation in a combined source. The resulting Boltzmann equation is 

Y\x) = -\x- n ~ 2 [Y 2 (x) - Y 2 q (x)] - S(x, cj>)Y(x)/x, (31) 

where 

S(x,4>)=<f>-Ct(9 + 2x)x- 2 * (32) 

and 

C, = 2. 9o 2 C 1 -^m^^(A 2 }o/M s 2 . (33) 

We are especially interested in the regime of temperatures mx 3> T; that is, as x —> oo. However, the asymptotic 
matching requires a knowledge of the solution for higher T as well. From (|32|) it is clear that for the case = 1/2 the 
x-dependence of the source and sink coincide for large x. Just as in ([5]), it is convenient to rewrite pip in the form 
of a differential equation with the nonlinear terms on the right side. We introduce the function 

g(x, 0) = x^ exp 

which is smooth at = 1/2, and the function 

f(x,cj>) = x- n - 2 /g(x,<f>). (35) 

We then define 



_9_ 

20 



2x 
1 - 20 



-20 



2C0 
1-20 



(34) 



Z(x,4>) = g(x,4>)Y(x), 

which is the analog of ((6|), and 

Z cq (x,<p) ee g(x,<p)Y cq (x) 

The Riccati equation satisfied by Z(x, 0) is 

dZ(x, 0) 



which is similar in structure to ([S]). The explicit form of f(x, 0) is 



= -A/(:z,0) \Z\xA) - ZLixM , (36) 



/(a;,0)=x-^ 2 ^exp 



9 2x \ _ u 2C 



^ { 20 1 



20 1 -20/ 1-20 



(37) 



The dominant asymptotic behavior of f(x, 0) as x — > oo changes according to the value of 0; f(x, 0) decays for 
> 1/2 and f(x, 0) increases exponentially for large x for < 1/2. Note that the phenomenologically relevant 
quantity is the Hubble-constant-free-relic abundance, ilh 2 = mN '/ 'pg, where pg is the critical density today and AT is 
the number density of the DM species. This is the quantity that is measured in experiments. For DM species X with 
mass mx it is given by [3j 

n x h 2 = m x Y(x)/x 3 . (38) 

The modification of fix can be compared to the standard (source-free) relic density by considering the phenomeno- 
logically interesting ratio 

n = lim „ lim (39) 

where J7^ uicc ^ frcc denotes the relic density of the DM species X in the standard cosmology case with constant dilaton 
and no space-time foam. We now systematically consider the behavior of the solution to the Boltzmann equation for 
various values of 0. 



A. The case of (f> near 1/2 

To investigate the behavior near <fi = 1/2 we let <j> = 1/2 — S and treat 5 as small. We write Z${x) = Z(x, 1/2 — 6). 
The Riccati equation satisfied by Z${x) in (131))) is 



Z' s (x) = -Xfs(x) [Zj{x) - Z 2 eq;S (x)] , 



where Z eq ^(x) = Z cq (x, 1/2 — S) for small S and 



with r)g = Ci/2— <s- Moreover, we have 



fs(x) « a; -™-5/2+«+2w exp (-gjig/x) 



Z eq ^{x)~Ax 2 - 5 ~ 2 ^e 



(40) 
(41) 
(42) 



As we did in Sec. [TT1 we argue that for large x in the post-freeze-out outer region, Zs(x) « Z| ost frcozo out (o;), 
where 



dx 



The solution to (gSJ) is 



r^post— freeze— out 



r _l _ , x - n -3/2+S + 2n s ' 

u <5 A n+3/2-(5-2?7 5 



and C5 is an integration constant to be determined. 

In the equilibrium outer region, following JS]), (J7J, and ([5]), we substitute 



(43) 



(44) 



Z s {x) ~J2\- k Z k , s ( 3 



k=0 



into (|4H|) . This leads to = Z eq $(x) and 



(45) 



Zl ^ (x) = -2^)i l0gZ - 5 - 
The value x = Xf, which characterises the freeze-out region, is determined by Zq g(x) = Z\ s(x), and we again obtain 
In the inner (freeze-out) region we introduce X as in (fT2| . The resulting equation for Z$(X) is 



1 d 



KdX -Z 5 (X) = -Ax- f 
The criterion of dominant balance requires that 



t»-5/2+*+2t)j 



Z 5 2 (X)-A 2 xy 25 ~^e- 2x f 



1 x — n— ^+S+2r)s 

— — AX j 



(46) 



(47) 



Following earlier arguments [see (1131) ]. in the inner region we have 

d 



dX 



Z 5 {X)^-Z 5 Z {X). 



The solution to this equation is 

Z 5 (X) = l/(X + V s ). 

where T>$ is a constant of integration. Matching (|44"|) with (|48|) gives 

1 A 



C & (n + 3/2-5- 2r) S ) 



n+3/2-5-2r]s 



(48) 



(49) 
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As in (fT9)) . the matching of the solutions in the equilibrium and freeze-out regions determines that 
The analog of @ is 



V s = Xx- n - 5/2+2ru+s . (50) 



1 _ n + 3/2 - 2?7^ - 5 + xj _ n _ 5 / 2+2%+(5 

C 5 ~ ?i + 3/2 - 2% -5^/ • l ° iJ 

Here, Xf = Xf ins . Consequently, for <f> = 1/2 — 6 and <5 small, the large-x asymptotic behavior of Y(x) is 

Y(x) ~ Cja;- 1 /^^. (52) 

In such a case and the freeze-out region are determined from (|5ip and from (jlOp . while the freeze-out point is 
given by (fTT|) . We note that the limit 5 — >• is smooth. 

B. The case of general (f> > with not near 1/2 

By the arguments given in Sec. HI] for large x in the post-freeze-out region the approximate solution to (|36|) is 

Z(x,4>) w ZP° st - frcozo -° ut (x,<?!)), where 

^_^post-free Z e-out( a . ) ^ = _^( a!j 0) ^post-freeze-out^ 2 _ (53) 



The solution to this equation is 



Z post-froczo-out/ = / 54 ) 



where is a positive constant. Equation (|54|) is valid for general cf> > 0. 
It is convenient to rewrite and (j3"T)l using the function 



We then have 



h{x,<t>) = {l-x 1 -^)/{l-2<i>). (55) 



(;(*, 0) - x* exp C ^V 2 ^) exp [2^h(x, </>)] 



and 



In the limit as x —> 



20' 

/(x, 0) = x- n - 2 -* exp f-^aT^ exp [-2&ft(*> <f>)} . 

( ' 0) l-2^T. for 0> 1/2, 

and so f{x, 4>) — > 0. Furthermore as x — > oo, Z cq (x, 4>) is negligible because in this limit 

f -oo for 0< 1/2, 
MX ' 0) ^l-2^T for ^> 1/2. 

1. The case < C<j> < 4> < 1/2 

Next, we consider the case for which £^/0 ~ 0(1) and C^x/ ^ 1- This case illustrates the competition between 
space-time foam and dilaton sources in their effect on the relic abundance. In this case 

g{x,4>) ~ a;* exp (j^J > 

f(x,<f>) ~ x-^-^exp^V (56) 
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for x in the freeze-out region and x Xf. The analog of ([S]) is similar except that A is replaced by Aexp ^— ■ 
The analog of CEU]) is 



x f ~ log 



(n + 1/2) log = log(2AA) - (n + 1/2) log (x f ) - ^. 



The previous analysis then implies that for large x we have the following asymptotic behavior for Y{x): 



where 



Y{x) 



x f = (n+l/2)W 



(n + 1 + <P)x n f +2+ ^ 
X(n + l + (f> + Xf) 



1+1/2 



(n + 1/2) 



(57) 



(58) 



(59) 



We denote this value of Xf by at/ i. The scaling (|58p is formally similar to the pure time-dependent dilaton case in 
(|23| , but the effects of the D-foam are incorporated only in the shifted value of the freeze-out point Xf in ({59} . 



The case > 1/2 

For </> 3> 1 /2 and a; > we have 

z = ~ Aexp U^lJ^TT+7 + c< (60) 

where C is a constant. To leading order the analog of (O for this case is independent of <f> and C</>, so is determined 
by (I10p . In the inner (boundary-layer) region we again write x = Xf + and iJ (X) = Z + kX). Hence, 

where (a;/) = -^ix^ 2 ^ 2n exp ^— 2 ^ 1 J ■ The principle of dominant balance then implies that 

1 \ -n-2-<4 / 2^0 

— = A x f exp 1 



20-1 



Hence, 4j| = — iJ 2 with the solution Z{X) = l/(x + T>), where D is a constant. 
Matching the equilibrium region to the boundary layer gives 

2Z cq (x f + kX) ps l/(x + T>). 

This implies that 



P = A. r - 2 -exp^j. (61) 
Matching the freeze-out- region solution to the post-freeze-out-region solution (|60p , we find that 

h-^{^i) x ^*( i+ ^h*y <62) 

Finally, we find that as x — > oo, 

Y{x) ~ i exp f - -1^-") ? — , . (63) 

A V 2^-17 ^-"-2-* + f x-™- 1 -* - x- n -^-A /(n +1 + 0) 

and x/ = Xf >ns . 



2& 
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3. The approach to <j> = 

The integrating factor in ([34]) is singular as <p — > + . However, the function g is only determined up to an x- 
independent factor. In order to study the limit <fr — > we consider a modified g(x, 4>) and an associated f(x, (f>), which 
we denote ])(x, <p) and f(x, 4>) respectively. These functions have the following form: 

g(x, q>) = x^ exp [2C,<f,hi{x, (j))] exp [^^(x, <p)} , (64) 

where 

l_ x l-2<l> X ~2<t>-1 2 

hi(x,4>)= l _ 2(j) , h 2 (x,<p) = — -, (65) 

and, as before, we have the relation 

f(x,cj ) ) = x~ n ~ 2 /g(x^). 

The limits x — > oo and <f> — ¥ do not commute (a feature that is common to other limits involving x). Parallel to 
the discussion of Ref. we take the limit — !• first. It is straightforward to show that for large x but C,x still small 
one obtains 

l/C = V + \x] n - 1+% /( n +l-9C). (66) 

By matching the equilibrium region to the freeze-out region we obtain 

V = Xxj n - 2+9C . (67) 

These formulas are similar to the the case of dilaton cosmology in the absence of space-time foam with the crucial 
difference that <f> is now replaced by — 9£. Finally, we obtain 

x % 

Y( - X) ~ C- 1 -Ax-«-i+9C/(n + l-9C)' (68) 

which indicates the role of D foam as a source of particle production in this case, in the sense that Y increases as 
x increases. Also, in this case Xf = x/ >ns . Notice that the behavior (1571) . which indicates an increase of the DM 
thermal relic abundance with decreasing temperature, is compatible with our earlier numerical investigations in Q. 
For x Xf (as in the current universe) the abundance ([68]) can be approximated by 

Y(x) ~ \~ 1 Xj +2 (x/x f f C (x > x f ), (69) 
which we use in the Sec. IIVI to discuss the phenomenology of these models. 



IV. PHENOMENOLOGICAL IMPLICATIONS 



As mentioned earlier, the phenomenologically relevant quantity that can be compared directly with experiments 
is the Hubble-constant-free relic abundance, £lh 2 = mN '/ 'pg, where p% is the current critical density and M is the 
number density of the DM species. For DM species X with mass mx this quantity is given by O [3]. The behavior 
of fix is then readily obtained for all cases studied in this work. 

The analysis in the previous sections indicates that time-dependent sources in our cosmological models lead to 
modified relic abundances for DM species, as compared to those computed within the standard cosmology. This mod- 
ification can be quantified by considering the ratio (1391) in which the numerator and denominator may involve different 
freeze-out temperatures. Since both expressions are known theoretically, the ratio (|39[) is computable explicitly for 
all cases studied above. 

Before proceeding with the phenomenology of the various sources discussed in this article, we make some generic 
remarks. If the sources are such that there is dilution of DM relic abundance relative to the prediction of standard 
cosmology, this can have important phenomenological implications for new physics, such as supersymmetry (SUSY) 
at colliders [1, S @ . In such a case there is a larger portion of the available parameter space of the SUSY model, 
which is compatible with the WMAP and other cosmological/astrophysical data 
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More room for supersymmetry implies heavier partners, which in turn may have interesting signatures at colliders, 
such as the Large Hadron Collider (LHC). If the relic density of the neutralino Xi> which is the dominant DM in 
SUGRA-like models, is diluted by a factor of about 1/10 in the presence of sources, then the final states expected 
at the LHC consist of Z-bosons, Higgs bosons and r-leptons. Such states are produced when one looks at the decay 
chains of the dominant SUSY production mechanism of squark q and gluino g pairs at the LHC: 

Q -> QX2 -> QTTx -> qrTxl, x° -> h°Xx, X° -> z Xi, 

where x°> 1S the next-to-lightest neutralino, and ft. is the Higgs particle. In Ref. a detailed analysis in the standard 
parameter space mi/2,77io (where mi/2 and mo are the gaugino and scalar masses) of mSUGRA models has been 
performed. In this analysis the parametric regions for the dominant decay patterns at the LHC: 

1. Higgs+jets+missing transverse energy, 

2. Z+jets+missing transverse energy, 

3. 2T+jets+missing transverse energy, 

have been predicted. Dilution factors of about 1/100 or more are compatible with the analysis in this paper for 
reasonable values of the parameters. Such dilutions may even push the parameter spaces of minimal supersymmetric 
models beyond the reach of the LHC (assuming standard-modcl-likc Higgs particle masses of about 125 GeV). For 
instance, in the constrained minimal supersymmetric standard model (CMSSM) with Higgs-mass range 123 - 128 
GeV and tan/3 of about 50, Lahanas and Spanos discussed the dilaton-induced dilution factor 6]. On including the 
effect of the dilution factor, they showed that the constraints placed on the parameter space of CMSSM, from the 
current ATLAS and CMS SUSY searches for DM, were not sufficient to exclude the model. 

We proceed to discuss the phenomenology of the cases discussed above by giving the corresponding values of the 
ratio (|3"§f today. We assume that the freeze-out points x f <ns in the absence of sources are about 30 , as expected in 
typical phenomenological models in which the DM is identified as a supersymmetric partner, such as a neutralino. 

The temperature of the universe, which is used in the definition of x today xq, is that of the cosmic microwave 
background (CMB) temperature T C mb ~ 2.35 x 10~ 13 GeV. Thus, for DM masses in the range m x ~ 300 GeV - 
1 TeV, we have 

x = tox/Tcmb Pa 10 15 - 10 16 . (70) 

Moreover, we assume that the source-free relic abundance Y ns (x) currently, which approaches a constant as x — > 00 
Q, as the boundary-layer analysis of Ref. [l[ confirms, is given by 

lim Y ™ " w .ui l V (71) 

Recall that the freeze-out point in the source-free case x / >ns indicates a range of values of x in the vicinity of (|10p 
with A Pa 0.000145 Q. For all but the case 0<( c j ) <(j)<gil/2 the freeze-out point Xf — Xf tns . However, as is 
evident from (f5"T|). even in the case < < <f> <C 1/2, the freeze-out point is shifted by an amount less than 9/2: 

a_0<si/2 ~ Xf_ ns — In the models we consider here Xf tTls pa 30, so such a shift is not significant. Thus, from now 
on we treat Xt pa Xf lV , s in all cases. This simplifies the arguments and allows an easy estimate of the ratio 1Z in (1391) . 

As a starting point, we take the case of a time-dependent dilaton source of the form ([3]) in the absence of D-foam; 
that is, £0 = 0. This case was discussed in Refs. [!, [f| and was revisited in Ref. [lj using asymptotic matching 
techniques. From (|2"3"|) and (|7ip the ratio (|3"§)) becomes (upon setting Xf ~ £/,ns) 

^ dilaton (* = ,0) ~ T'lT n (72) 

n+ 1 n + l + <f> + Xf, na 

with Xq given in ((70)) . 

From (|5"8")) we then notice that (|T2"|) also applies to the case of nontrivial D-foam but with <(</>< ^ 1/2- For 
Xf^ sn about 30 and for s-wave scattering (71 = 0) the approximate thermal DM relic dilution factor (|72|) is determined 
by (x/^s/xo)* ~ 10~ 16 "^ for DM masses m x in the range 0.3 , 1 TeV. Thus, to obtain a dilution factor of order 1/10, 
which is relevant for LHC phenomenology, we need values of </> near 1/16, which is small compared with 1/2 and 
which is consistent. However, the case of phenomenologically significant dilution requires that Cip x f 1 an d thus 
(cf, <C <fi. For the pure dilation case, in the absence of D foam, one may have larger values of <j> that lead to acceptable 
phenomenology; for instance, a dilution of about 10 -2 can be obtained with <fi pa 1/8. 
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On the other hand, in the case where the space-time defect (D-foam defect) dominates the time-dependent dilaton 
effect, that is when the strength of the foam fluctuations is such that £ 3> (j) ~* 0, we have an enhancement of the 
DM relic abundances rather than a dilution as the temperature decreases. This becomes clear from (|68| and ()69|) . 
In such cases there is less room for supersymmetry available in the relevant parameter space as compared with the 
source-free case after cosmological (WMAP) constraints [7J are taken into account. 

The enhancement factor scales like 

1l(x = x Q )~ n+ n l + l Xf (x /x f ) 9i . (73) 

For s-wave scattering and with xq given by (f70|) this implies that 7Z ~ (lO)' 136 154 - >l \ Such models lead to more 
severe constraints on the available supersymmetry parameter space if the enhancement is observable. 

Therefore, for these models to be phenomenologically viable today, this requires 1Z to be 0(1) within experimental 
error, so that the increase compared to the source- free (standard) case is not appreciable. This requires that C < 10~ 3 , 
so that the error in calculating abundances would match the per mil level of the current errors in experimental 
astrophysical measurements Q ■ Because of (|3"3"j) , this implies that 

C = 2x (g m x /M s ) 2 (A 2 ) < 1(T 3 . (74) 

For the range of xo in ([70)1 this is satisfied for heavy D particles of masses M s /go > (lO 11 — 10 13 ) /(A 2 )oGeV, and 
m x in the range 0.3 - 1 TeV and (A 2 ) < 1. 

Next, we discuss the case of a time-dependent dilaton with <fi 3> 1/2 in (1531 . Since we always assume weak 
foam fluctuations, < 1) this case is also dilaton dominated, and hence we expect a significant dilution of the 
relic abundance. Indeed, because of (|63[) and (|71[) . in this case we obtain the following expression in the limit 
x — >• xq 3> Xf = Xf ns for the ratio 1Z in Eq. (|39[) : 

* ~ •* ( - 5&) (Hrr) C 11 ^) <»*>~r> (.//•)' . w 

where we assume that the freeze-out points Xf ~ Xf^ ns are about 30 or more. Thus, we have significant dilution of 
the DM relic densities at late epochs of the universe. For instance, in the present era and for DM masses in the range 

m x « 300 GeV, the ratio is K « cxp (-j^^A x (« 200) x lO" 15 ^ for s- or p-wave scattering (n = 1,2). Thus, we 

see that for £ -c 1, which is natural in the case of D foam with D particles whose masses are higher than a TeV 0], 
the main factor that drives the dilution is the value of the dilaton parameter <j>. 

In the case <f> « 5, for instance, the dilution factor is already enormous (it is of order 10 -75 ), so in such models 
practically all DM today will have disappeared. This may rule these models out phenomenologically, although the 
situation with DM and its nature is currently unclear, as there is no concrete evidence for it apart from the galactic 
motion. For this reason alternative theories with no dark matter but modified gravity at galactic scales have been 
considered extensively in the literature. We do not consider them here, since in our opinion the evidence against them, 
especially from galactic lensing measurements, is significant. Thus, all we can say is that this type of supersymmetric 
DM (satisfying Xf ~ 30) would be diluted in this model, and it would be practically absent today. Other types of 
DM that would not couple to the dilaton, might survive. 

Next we discuss the cases for which <fi is near 1/2. Now, using (f5Tj) and ([52")) , we obtain for x — >■ Xq 3> Xf ~ Xf tns : 

n ~ "t / 3 2 /2 ~ 2Cl/2 -f,~ g n + 1 l; / '" s (W*) 1 ^-'"' ■ (76) 

n + 3/2 - 2Ci/2-<5 + Xf,™ n + l 

We observe that for 6 — > and C1/2 < = 1/2, the main dilution comes from the dilaton effects and scales 
with x as (xf^ns/x) 1 / 2 ^^ 1 ' 2 . Thus the dilution due to the dilaton is compensated by foam fluctuation effects, so 
f° r C1/2 = 0/2 = 1/4 < 4> = 1/2 there is no appreciable dilaton-driven dilution, and the ratio ([75)1 tends to one 
(i?Ci /2 =i/4 ~ 1) for an y n > 0. 
We observe from ((33]) that d/ 2 is independent of xq: 

Ci/a = 2. g 2 (A 2 ) ™ 2 ,/M s 2 (77) 

and the condition £ = 1/4 implies that (A 2 )o = M 2 / (8<7o TUx) < lj where the inequality on the right side ensures 
naturalness in the fluctuations of a weak foam, which we have assumed throughout. The latter condition necessitates 
mx ~ M s /g s . We stress that in this case the result for the relic abundance today turns out to be equal to the 
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standard source- free case independent of the actual freeze-out point, and hence in principle mx is only constrained 
to be of the same order of magnitude as the D-particle mass M s /g . 

Finally, we mention that one may consider a 6 > 1/8 to produce dilution of order 1Z < 10~ 2 in the relic abun- 
dance (|76[) . thereby opening the possibility of pushing this class of supersymmetric models out of the reach of the 
LHC, according to the analysis in Ref. [6(. However, in this case, (|33l) implies that the condition 1/4 = £5 ~ 
2g$XQ S (A 2 ) m x /M2 for x in the region ([70]) can be satisfied for g$m x /M% ~ (1/8)10 15<5 (A 2 ) . To ensure that 
mx < M s /go this would imply naturally small fluctuations in the foam (A 2 )o ~ 10~ 15<5 with 5 > 1/8. 

The above predictions are quite generic and hence they are largely independent of the details of the underlying 
microscopic model. Nevertheless, the cosmology of the models, in particular the precise dependence of the dilaton 
on the cosmic time at various eras of the universe, is an open issue. The lack of detailed microscopic models that 
would determine the form of the dilaton potential and provide rigorous information on the region of validity of the 
dilaton cosmological solution ([3]) complicates matters. Nevertheless, one may perform phenomenological searches 
on the compatibility of such solutions at various epochs of the universe. For the DM searches mentioned above, 
all one needs is the dominance of the time-dependent dilaton at early epochs of the universe before the big-bang 
nucleosynthesis. Nevertheless, a dilaton of the form ^ can be compatible (notably at the same level as the ACDM 
model) with cosmological data even at low redshifts of order z = 0(1), where large scale structure in the universe 
(galaxies and clusters of them) is formed, as demonstrated recently in Ref. [T3|. On the other hand, D-foam dominance 
at late eras (such as the end of radiation or matter-dominated era [3]) has been argued to play a role in galactic 
growth itself. Thus, considering models with combined dilaton and D-foam sources, as in the current article, may 
be desirable from the point of view of constructing realistic cosmologies in such frameworks. However, the rate of 
galactic growth is in principle capable of discriminating the various models ([3]) corresponding to different values of <p 
when more data become available in the near future. In all such theories, of course, an important requirement is that 
the big-bang-nucleosynthesis conditions at MeV temperatures are not disturbed. 
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APPENDIX: Thermodynamic properties of a universe in the presence of sources 

The purpose of this appendix is to demonstrate that it is possible to define an appropriate entropy density (scaling 
with temperature as T 3 ), even in the presence of nontrivial backgrounds, such as a time-dependent dilaton and/or 
space-time foam. This allows the entropy density to be used in this paper as a fiducial quantity in the definition of 
the thermal relic abundance Y(x). 

In the presence of such nontrivial backgrounds the continuity equations of cosmic fluids corresponding to matter 
and radiation are modified relative to standard Friedmann-Robertson- Walker (FRW) cosmology. These modifications 
could affect the thermodynamic properties of the universe, such as the relation between the scale factor and the 
temperature T (the cooling law). It is the relativistic degrees of freedom that dominate the entropy and the cooling 
law. In the case of a FRW universe, the continuity equation is the conservation of the stress-energy tensor V M T M „ = 0, 
which is compatible with Einstein's equations of general relativity and admits a thermodynamic interpretation. This 
equation can be manipulated to appear as the first law of thermodynamics for the total internal energy pV in a 
co-moving volume V ~ a 3 and pressure p: 

d(pV) + pdV = 0, V ~ a 3 . (78) 

This interpretation in terms of the first law is consistent with an adiabatic expansion at temperature T. A constant 
entropy function S(T, V), analytic in T and V, can be constructed: 

S = vB+£. (79) 



The construction involves the application of the thermodynamic Maxwell relations to cast the right side of (l78l) into 
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the form TdS; that is 2 , 

= TdS = d{p V) + pdV. (80) 

The dominance of relativistic degrees of freedom in the entropy S is therefore consistent with the constancy of S and 
the cooling law a ~ 1/T. (Recall that p = 3p ~ T 4 for radiation.) From (|79[) and (|80|) it is then straightforward to 
see that the entropy density s = S/V scales with temperature T as T 3 . 

For cosmologies with nontrivial time-dependent dilaton and/or space-time D-particle foam backgrounds we also 
construct entropy functions that are constant during the evolution of the universe. 



(i) Dilaton Cosmology 

For a FRW cosmology in four space-time dimensions the presence of a rolling dilaton leads to a modification of the 
continuity equation for the energy density p and pressure p 



p + 3H(p + p)-<i>(p-3p) = Q, (81) 

where the dot denotes derivative with respect to the cosmic time t. (Here p and p here denote the total energy density 
of the fluid.) We note that: 

• From (|81[) the dilaton source terms do not play a role for radiation; one obtains the standard scaling of p ~ a~ 4 
in the radiation dominated era of the universe. 

• For dust, p = 0. Also, for DM with mass mx, Px = mxnx, where nx is the number density; the source- 
independent part of (fBTj) yields the collisionless Boltzmann equation for thermal relic abundance. The dilaton- 
dependent term is a classical source term <i>n. 

The nonlinear part of the Boltzmann equation comes from two-body annihilations of DM particles. On assuming 
the functional dependence p = p(a) and a dilaton source of the form ([3]), we obtain from (IcTTj) 



d(pV)+pdV- |3> |(p-3p)dV73 = 0. (82) 

Here we have taken into account that $ = — |$o|-ff, where H = a/a is the Hubble parameter, V ~ a 3 is the co-moving 
volume, and pa 3 is the total (internal) energy in that volume. 

We thus observe that the presence of a rolling dilaton affects the standard thermodynamic properties of the FRW 
universe. Our aim is to ascertain whether the total entropy in the co-moving volume V remains constant in time 
after the inclusion of the dilaton source ©. A naive application of the first law of thermodynamics would identify 
d(pV) +pdV with TdS, where T is the temperature, and S is the total entropy in the volume V. It would seem that 
a dilaton source leads to the nonconservation of entropy. However, this is incorrect. To show this, we first replace the 
zero of the right side of Q82j) by TdS, where S is the quantity that represents the entropy; the entropy is assumed to 
depend on T and V, so S = S(T, V). From fl52J| we find that 

dS = ^d(pV) + L dV -M{p- 3p)dV = V^dT + i [(1 + |$ |)P + (1 - |*o|/3)p] dV. (83) 

As in the case of conventional cosmology, it has been assumed that the cooling law of 1)2 7p holds, and that p = p(T) 
and p = p(T), since both p and p depend upon the scale factor a, which is a function of temperature, a — a(T). The 
function S(T, V) is assumed to be a differentiable function of T and V. This implies the condition: 



d 2 S d 2 S 



dTdV dVdT 

From and (EH we then obtain 



(84) 



^[d + |*0|)P + (1 - |*o|)/3)p] = { ^^% - (85) 



2 Strictly speaking Eq. 1 180 II should be modified by a term involving the chemical potential ^t; hence TdS should be replaced by TdS + fj,dN . 
As is standard in cosmology fi is ignored because fi/T is much smaller than one [3j, which is consistent with the dominance of the 
relativistic degrees of freedom in the entropy. 
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The expression ([83]) for dS can be rewritten as 



dS = Ld[(p + p)V}-^dT-^d[(p-3p)V] 



\$o\V dp dT |$ |V dp dT 
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In the last equality on the right side we have used (|85|). From ([86]) we conclude that the quantity 

s(t,v) = [ P (i ~ |$ |/3) + P (i + |* |)]vyr 



3T dT 

(86) 

(87) 



is constant upon using the classical equations of motion [or equivalently, the continuity equation (]8ip for the case of 
dilaton cosmology ([5])]. S may be identified with the total entropy in the co- moving volume V. The corresponding 
entropy density s is then: 

* = ^[p(l-|*o|/3)+i>(l + l*o|)], 



which, in view of (f87| . scales with the size of the universe as a~ 3 = (T/Co) 3 , upon assuming (f27|) . We stress that the 
energy density p and pressure p in the above formulas pertain to the total degrees of freedom of the fluid, including 
the relativistic ones. It is the latter, for which the dilaton source effects are irrelevant [see (|8"Tj) ]. that provide the 
dominant contributions to the entropy; otherwise, the entropy would not remain constant. Indeed, in the case of DM 
dust, p = 0, the entropy density is s = p(l — |<J> |/3)/T, which does not leave the entropy function ([57)) constant. 
This is satisfied only for relativistic degrees of freedom that have an energy density scaling like p ~ T 4 with the 
temperature T. 

We have the following relation between Y and the number density nx of the DM species X: 



Y = n x T~ 



m\Yx 3 , 



(89) 



as in the standard cosmology case. The energy density px of the DM relic satisfies px = ^ix n x- The current relic 
abundance: 



<-x> 



r 3 ' 

Po x 



(90) 



occurs for x — xq — mx /7b, with To the current (CMB) temperature of the universe and pg the current critical density. 
Since the latter is proportional to h 2 , the above expression is independent of the value of the Hubble-constant. 

In practice, io > 1. Hence, the asymptotic regime Y(x) with x — > oo is relevant. In the current literature one 
usually replaces Y by F^; that is, 



n x h 2 



m X 



Pa 



r.3 - 



(91) 



For standard cosmology in ([6"2[) lim x _ J . 00 Y(x) — constant. This constant value of the freeze-out is identified with the 
current relic abundance of the weakly-interacting massive particle (WIMP) flxh 2 ~ l/(av). 

We remark that the scaled Hubble-constant-independent relic abundance of the DM species X behaves as 



n x h 2 



m^x Yrx 



(x 



Also, for the case of standard cosmology Foo = constant [see ([62]) ]. 

fl x h 2 — x~ 3 ->• (x -> oo). 
For dilaton cosmology [see ((3])] one has a modified law 



Cl x h 



x 



-3-|* | 



-4 (x — > oo) 



(92) 



(93) 



Finally we remark that for the dilaton case $o > 0, the string coupling would increase for large times, and the 
theory would become strongly coupled and thus intractable. Nonperturbative string corrections would need to be 
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incorporated. Nevertheless, the formal solution for Y(x) behaves asymptotically as Y{x) ~ x'* '. This would still 
imply an asymptotically vanishing relic abundance provided that $o < 3 because we have Qxh 2 ~ x~ 3+ l*°l — > as 
x — > oo. 

(ii) Stochastic D-Particle-Foam Cosmology 

It is known that in the background of D-particle space-time foam (for constant dilatons), the Boltzmann equation 
for the number density nx of the DM species X assumes the form [9J] 

^-n x + 3Hn x = r D _ foam (i)n x + C[n], (94) 
at 

where Tn_{ oam (t) = 2Ha 4 mx < A 2 > jfzT (9 + ^ m x/T). The notation and conventions here are those of where 
C[n] = (o~v) [(n°x) 2 — (nx) 2 ] is the standard nonlinear interaction term and n ^ is the thermal equilibrium number 
density of X. As discussed in Q and reviewed in the text [see Eq. ([^6]) ]. the recoil fluctuations of the D-foam <C A 2 ^> 
averaged over populations of D-particle defects, have the scaling <C A 2 (A 2 )oa~ 3 . 

To this end, we use the cooling law (l27t and ignore the nonlinear interaction term C[n). From (|94|) and (|26|) . for 
the regime of temperatures rax 3> T, the energy density px = mxnx then satisfies the continuity equation 

■^px + 'ZHpx = f D ^{ oam (t)Hpx, (95) 

~ 2 2 ~ 

where rrj_f oam (i) — 4Co-^-^ L (A 2 ) . For weak foam effects we have rrj_f oam < 1 in the range of temperatures we 

are interested in; that is from the early universe until today (T > Co). Equivalently, for an expanding universe where 
a > we have 

d(p X V) - pfD-foa.ndV/3 = 0. (96) 

Equation (|96|) implies that the thermodynamic interpretation of heavy DM dust in the foam background is that of 
a gas with an adiabatic expansion of its volume. During the expansion entropy is constant, and the effective pressure 
Pcff-x of the gas is negative (indicating cosmological instabilities): 

Pctt-x = -pr D _ foam /3. (97) 

Note that p c s~x has a nontrivial dependence on the temperature. 

From the cooling law (1271) . we may write (see [SS]) T = ja, where 7 is a constant much less than one. Hence, the 
scaling of the dust energy density, due to its interaction with the D-foam, is easily obtained from (|96[) to be (in units 
of 00): 

px - a- 3 K da ~ TV i, ^- 1 ) (98) 

To find the entropy function that remains constant it is essential to consider the total energy density p, including 
relativistic degrees of freedom, and not only px- In a similar spirit to the dilaton case, the relativistic degrees of 
freedom are insensitive to the heavy D-foam source effects; in this sense they satisfy an equation of the form (|78j) by 
themselves with equation of state p = p/3, which can be added to (|96|) to give the equation 

d((p rad + PX )V) + (p rad + PcS - X )dV = d(pV) + PcS dV = 0. (99) 

Equation is the analog of the continuity equation in the case of D foam. We stress that in (|99p p and p e s refer 
to the total energy density and pressure, including relativistic degrees of freedom and D-foam background effects. 

Taking into account that p is a function of T, we can formally replace the right side of (|99| by TdS to determine 
the (constant) entropy function S (ignoring chemical potential terms, a valid assumption for weak D-foam); only at 
the very end of the computation will we set dS to zero. We then have 

dS = Y.^. dT+ E+m dV , (100) 

T dT T v ' 

S is considered to be a smooth function of T and V , which are treated as independent variables. From the requirement 
we deduce the condition 

_ P+PeS jjPeg , . 

T 2 T dT { ' 
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We then see immediately from (jlOOl) and (|101[) that 



dS 



P+Peft 
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Vdp, 



'err 



T dT 



dT 



PcS 



VdT 



PcS 
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V 



which upon setting dS = implies the constancy of the effective entropy function in the co- moving volume V: 

S = S cS = P+ ^° S V = constant. (102) 

Note that we have used ([98| and the cooling law (|27|) : p and p e ff refer to the total energy density and pressure 
including rclativistic components, which is essential for consistency. As in the previous cases, the relativistic degrees 
of freedom dominate the entropy. The entropy density s associated with S is given by an expression similar in form 
to that in standard cosmology: 



Seff D-foam = P + Pcs/T 



(103) 



and scales with the temperature as T 3 . Hence, s can be treated as a fiducial quantity to define Y(x) just as in the 
dilaton cosmology case (i) above. 

Notice also that for the case of dust in dilaton cosmology, the effective entropy function (|88|) is reproduced upon 
replacing the source To-foam by the corresponding source of the running dilaton ([3]) cosmology r runn ing dil = — l^ol- 
(With our definitions we have r running du = $ = -|$ |ff = f mnning dii H.) 

In the paper we considered the combined source case, where the foam appears together with a nontrivial running 
dilaton of the form ([3]). The string coupling g s — e* exhibits a nontrivial scaling with the scale factor and also with 
temperature. The combined source is taken to be the algebraic sum of the respective two source terms; that is, 



total 



- total 



H = 



l-2|$o| 



H 



(104) 



in the asymptotic region rax 3> T of interest, where we have assumed the cooling law (127 
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